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verse. We consider the stochastic evolution of density parameter in the early universe 
after the inflationary phase qualitatively, under the assumption of fluctuating w factor 
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the universe depends on the energy density, from the coupled Friedmann equations 
we calculated the two variable probability distribution function assuming a flat space 
geometry. 



Key words. Cosmology; stochastic equations; Fokker-Planck equation. 



PACS Nos 98.80; 02.50 



f e-mail: sivakumarc @ cusat. ac. in 

| Permanent Address: Department of Physics, St. Thomas College, 
Kozhencherri, Kerala, India 689641 e-mail: money @ stthom. ernet. in 



1 



1. Introduction 

The simplest model of the universe [1], called the Friedmann model or hot big bang 
model, is based on the assumption that the matter distribution in the universe is ho- 
mogeneous and isotropic on very large scales (cosmological principle) . This assumption 
simplifies the problem of solving Einstein's equations. The assumption of homogeneity 
leads to the prediction of deterministic Hubble's law and one expects a scatter free 
redshift-magnitude diagram for galaxies and other extra galactic objects. For nearer 
galaxies, Hubble's law is justified by a relatively scatter-free redshift-magnitude rela- 
tion. This corresponds to a deterministic evolution of the universe. However, for high 
redshift quasars (in the Hewitt-Burbidge catalogue [2,3]) and for the latest Type la 
supernovae data [4,5], the redshift-magnitude diagram is a scatter diagram, i.e., they 
are in contrast with the deterministic Hubble's law. If Hubble's law is valid for all 
extra-galactic objects, then the evolution of the universe (or expansion rate) is non- 
deterministic in the early epoch, as indicated by the scatter diagram. We observe 
that the scatter increases as we probe into more and more distant epochs. The small 
scatter for galaxies at low redshifts is explained as due to peculiar velocities and the 
conventional explanation for the peculiar velocities is that they are induced by the 
observed density perturbations. This may be adequate to account for the observed 
peculiar velocities of objects, in the range of 100 km/s - 400 km/s. Thus this may 
lead to a very small scatter at low redshifts. But since the amplitude of density per- 
turbations in the early universe was very low, the large scatter at very high redshifts 
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remains unexplained and it is desirable to look for some alternative mechanism. Some 
authors speculate [3,6] that the quasars are not at their cosmological distances and 
have proposed some non-cosmological contributions to z as a possible explanation for 
the scatter in the quasar data. Since the distance measurements are extremely difficult 
for quasars, due to the difficulty in identifying standard candles, the scatter may be 
caused by the variation of intrinsic luminosities of quasars of same z. However, the 
data is more accurate for supernovae and the scatter in its redshift-magnitude diagram 
is not due to either peculiar velocities or variation of intrinsic luminosities. In [7], 
we proposed that a fluctuating equation of state or a fluctuating mean w factor in 
the equation of state led to a non-deterministic or stochastic Hubble parameter and 
argued that such a fluctuating expansion rate in the early universe might have led to 
a randomness in the recession velocities of objects, in addition to peculiar velocities 
and will produce a scatter in the redshift-magnitude diagram in those epochs. Here 
we develop a more general description of the stochastic dynamics of the early universe, 
and discuss the non-deterministic character of total density of the universe, as well as 
the scale factor for the universe. Such a stochastic approach is necessary, when the 
mean w factor in the equation of state of the cosmic fluid is a fluctuating quantity. 

2. Observational issues and assumptions 

In this section, we state the conditions under which a stochastic equation of state 
(or a fluctuating mean w factor) emerges. In standard cosmology, the cosmological 
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fluid is in fact not unicomponent, instead matter and radiation (with equation of state 
p m = and p r = (|) p r , respectively) in disequilibrium coexist in many 'elementary 
subvolumes' of the universe [8]. Some recent measurements on the age of the universe, 
Hubble parameter, deceleration parameter, gravitational lensing etc., point to the need 
of extending the standard model by including some new energy density (missing energy) 
in the present universe, in addition to the usual relativistic/non-relativistic energy 
density Recent observations using Type la supernovae [4,5] as standard candles seem 
to indicate that the universe may be accelerating, driven by a positive cosmological 
constant (or a vacuum energy density p v ) with the equation of state p v = —p v . Some 
authors [9,10] introduced Quintessence or Q-component (for exa: scalar fields rolling 
down a potential with negative pressure or cosmic strings etc.) with equation of state 
Pq = wq pq , — 1 < wq < 0, in addition to matter density. The case for a positive 
cosmological constant has been considered very often in the literature [11-23]. All these 
models favour a flat universe with Q = Q m + Q v + ■ ■ ■ = 1, where Q is the ratio of 
energy density to critical density. Thus our universe may be approximated by a perfect 
fluid having many components, each with equation of state pi = WiPi, — 1 < Wi < +1 , 
i — 1, 2 • ••. If we denote the total energy density due to all such components as p, then 

P = Pm + Pr + Pv H , (1) 

where p m , p r , p v etc. are the average densities of matter, radiation, vacuum energy etc. 
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In a similar way the total pressure p can be written as 

P = Pm+Pr+Pv^ (2) 

In general, p — J2Pi, where pj's represent energy densities of various components. From 
the energy-momentum conservation law (here it is assumed that only the total energy 
density is conserved), we have 

p = -3-(p + p) = -3-p(l + w), (3) 
a a 

where p = J^Pi is the total pressure, and the ratio w = p/p should lie between -1 and 
+ 1. Splitting p and p into individual components, the above equation becomes 

Pi + P2 + Ps + • • - = -3-[pi(l + Wi) + p 2 (l + w 2 ) + p 3 (l + w 3 ) + ••]. (4) 

(X 

From Eqs. (3) and (4), we get 

W (t) = - 1/3l+ /^f ' ] -l= 1Pl(1 + Wl)+P2(1 + W2) + " ,,] -l. (5) 
v ' 3(a/a)p p y ' 

Here we assumed that the total energy density p is conserved and not the parts corre- 
sponding to pi,p2 • •• separately, and hence there can be creation of one component at 
the expense of other components and since a/ a and p are large in the earlier epochs, at 
least some of the /Vs will be significantly large in this period. Since recent observations 
indicate the existence of vacuum energy even in the present universe, one can expect 
that particle creation continued for a fairly long period in the early universe [24]. If we 



have a many component fluid, then the Einstein equations, along with the equations 
of state of individual components, are insufficient to determine the creation rates of 
each component. In [25], Weinberg discusses the consequences of the presence of a 
cosmological constant (vacuum energy) in the energy density. He discusses some phe- 
nominological proposals made by some authors, of the energy transfer between vacuum 
and matter or vacuum and radiation, in such a way that either p v j p m or p v j p r remains 
constant, respectively. He also considers the possibility of creation of radiation from 
vacuum energy, keeping p v /p m fixed. However, in a general case as in Eq. (5) above, 
one cannot expect any kind of creation in the universe to be a smooth process, since 
they can be sporadic events occurring in different locations and times, like those occur- 
ring in galactic nuclei. Here, as in the case of other stochastic processes like Brownian 
motion, a complete solution of the macroscopic system (universe) would consist in solv- 
ing all the microscopic equations describing the creation processes, but such a rigorous 
derivation will be very complicated or even impossible. In this context, a stochastic 
approach is more reasonable, in which we consider the creation rates to be fluctuating, 
leading to fluctuations in the ratios Pi/p. As it clear from Eq. (5), this in turn, will 
lead to a fluctuating w factor, which is the key assumption made in this paper. Such a 
fluctuating quantity modifies the dynamics of the early universe, where the evolution 
of the cosmological parameters, like the total energy density (p) and scale factor for 
the universe (a) becomes stochastic or non-deterministic. The fluctuations in the ratio 
Pi/ p, that we are taking in to account here are classical, i.e., our stochastic model is a 
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modification to the classical Friedmann model of the early universe, when fluctuations 
are significant. In [22], Fang et al discuss a stochastic approach to early universe (before 
recombination epoch), when cosmic fluid consisting of primeval plasma and radiation, 
is not perfect, but have dissipations due to differences in the adiabatic cooling rates of 
the components of the fluid and the possible energy transfer between them. Physically 
motivated interaction models are also proposed in the literature [26,27], which lead to 
energy transfer between various components. However, once we probe into still earlier 
epochs (stages of inflation etc.), quantum fluctuations become very important. Many 
authors discuss the need for a stochastic approach to inflation [28,29,30], when the 
quantum fluctuations of the scalar field are significant, and try to get a probability 
distribution function for the scalar field after solving the quantum Langevin equation 
(or FPE) describing the evolution of the scalar field. However we adopt the stochastic 
approach in the classical regime, where fluctuations in the creation rates and also in 
the possible energy transfer between different components of the cosmic fluid lead to 
a stochastic equation of state. This causes a non-deterministic (stochastic) expansion 
rate for the universe, described by a set of stochastic differential equations, instead 
of the deterministic Friedmann equations and we evaluate the probability distribution 
function of the cosmological parameters. In section 3 we discuss the stochastic charac- 
ter of the density parameter of the universe on the basis of Fokker-Planck formalism 
[31] . In section 4 we analyse the effect of such phenomena on the expansion factor for 
the universe and a two variable distribution function is derived. 
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3. Stochastic evolution of density parameter 

Suppose the universe is approximated by a many component fluid in the early epochs, 
with a fluctuating w term in the equation of state. Now we write the evolution equa- 
tion for the total density in the early universe (assuming that, total energy density is 
conserved), immediately after inflation, when curvature factor appearing in the field 
equation is negligible, so that the background is approximately flat. 



p = -?>-[! +w{t)]p. (6) 
a 

Overdots denote time derivatives. Using Friedmann equations we have 



p = -V247rG[l + w(t)}p 3/2 . (7) 

Above equation is a stochastic differential equation of the Langevin type. Since w is a 
fluctuating 'force' term, p is a stochastic variable or its evolution is non-deterministic. 
The random behaviour of p in the early universe is due to fluctuations in the factor w 
alone. If fluctuations are zero we are back to the deterministic standard model. We 
apply stochastic methods [31] for the analysis of the above equation and the probability 
distribution function is calculated using Fokker-Planck formalism. By making use of 
the transformation 

° = (GkGp)W> (8) 

eq. (7) gets modified into 

a=l + w(t), (9) 
8 



which is a non-deterministic, stochastic first order differential equation (Here a oc t for 
a pure deterministic case). To solve eq. (9) we use certain simplifying assumptions 
that the fluctuating factor w is Gaussian (^-correlated, with mean zero. Though these 
assumptions are taken for the sake of simplicity, we expect that they are reasonable 
when compared to the time scales involved. 
If we have a general Langevin type equation of the form 

y = h(y,t)+g(y,t)T(t), (10) 

where T(t) is a fluctuating quantity with zero mean and Gaussian (^-correlated, then 
the corresponding Fokker-Planck equation (FPE) describing the time evolution of the 
probability distribution function W(y,t) can be written as 

where D^ n \y) are the Kramers Moyal expansion coefficients given by 
D {n) (y,t) = 1 [lim - ([y(t + r) - x] n ) 



(12) 



y(t)=x 

Here y(t + r) (r > 0) is a solution of eq. (10) which at time t has the sharp value 
y(t) = x. Under the assumption of 5 -correlation and zero mean of T(t), all coefficients 
vanishes for n > 3 and retain only the coefficients and D^ 2 \ called drift and 
diffusion coefficients respectively. Following the standard procedure [31], eq. (9) leads 
to the FPE 

dW(a,t) dW d 2 W 
9 



Here we used drift coefficient = 1 and diffusion coefficient D^> = D (D is a 
constant with dimension of time, which is introduced for the purpose of generality). 
In order to obtain non-stationary solutions of eq. (13) we use a separation ansatz for 
W(a,t) 

W(a,t) = (f)(a)e- xt . (14) 
Substituting this into eq. (13) and solving for 0(<r) we get 



<j){a) = Aexp[^ + ika], (15) 

where 



Thus we see that for A < 1/4.D, k 2 is negative and the solution is exponentially diverg- 
ing, which is not a physically reasonable solution. Hence we conclude that A > 1/4.D, 
so that k is real. We write the most general solution as 

W(a,t) = J2c n <t>n(v)e- x " t , (17) 

n 

where c n can be real or complex but W (<r, t) is always real. For a continuous parameter 
k, from eqs. (14) and ( 15) the general solution or the distribution function is given by 

/ + CO fj f 

exp[— + ika-k 2 Dt-—]dk. (18) 
-oo ZU ^±LJ 

We choose A = l/2n for normalization purpose. On evaluating the integral we find 
the distribution function 



10 



ADt 



(19) 



which is in Gaussian form. The expectation value of the stochastic variable a is (a) = t 
and corresponds to the deterministic solution of eq. (9). The width of the Gaussian is 
found from the variance, v = ((a — (cr)) 2 ^ = 2Dt. Once W(a, t) is known it is straight 
forward to write the distribution function W(p, t) as 

1 



W (p, t) 



: CXp 



1 -ty/EnUp 
2AnGDpt 



(20) 



^/96n 2 DGpH 

We can also find the transition probability for the stochastic variable to change from 
an initial state (<r',t') to a final state (<x, i) as 



P{a,t\ a', t') 



sJAkD (t - V) 



exp 



[(a-a')-(t-t')f 
AD (t - 1') 



(21) 



with the initial value 



P(a,t | a',t) =5 (a -a'), 



(22) 



indicating Markovian nature of the random variable a. In terms of p eq. (21) becomes 



P(p,t\p',t') 



^AttD (t - V) 



exp 



pt-^-V^GWit-t') 



(23) 



247rGpp'D (t - f) 

This represent the probability for the energy density to change from an initial value p 
to a final value p' during a time interval (t — t') in the early epochs. This characterises 
the stochastic behaviour of density evolution in the early universe. 



4. Scale factor as a stochastic variable 



Under the assumption that the factor w is fluctuating during the early epochs, the 
evolution of the scale factor also becomes non-deterministic, since the time evolution 
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of a (t) is determined by the total density. So we have a system of coupled stochastic 
differential equations derived from Friedmann equations 



I8itG 



-P, 



(24) 



and 



p = -V24nG [1 + w{t)]p 3/2 . (25) 

Here we are considering the dynamics of the universe immediately after inflation, so 
that the background can be treated approximately flat. With the transformation de- 
fined in eq. (8), the above system of equations reduce to 

4 - !■ < 26 > 

and 

a = l + w(t). (27) 

Following the standard procedure [31] we have the drift coefficients = = 1 

and the diffusion coefficient D£> = D is assumed to be a constant. It shall be noted 
that this diffusion term arises due to fluctuations in w alone. The two variable FPE 
for the distribution function W (a, a, t) can be written as 



dW 2 



dt 3a 



W + a— 
oa 



We can solve the FPE by first assuming the ansatz 

W (a, a, t) = U(a)V(a)e- xt , (29) 
12 



and substituting into eq. (28). We obtain 



<j d 2 V adV 2 



a dU 
U~da 



+ 1 



(30) 



Each side in this equation can be equated to a constant m. When m = 



U(a) oc -, 
a 



and 



a 



V(a) oc exp[— + ika], 



(31) 



(32) 



with k given by eq. (16). A physically reasonable solution exists for A > 1/4D, which 



is 



B 



a 



2D 



+ ika — Xt 



(33) 



W (a, a, t) — — exp 
a 

Here B is a normalization constant, chosen to be l/2n. One point to be noted is 
the most general solution to eq. (30) when m ^ 0, is a series solution owing to the 
singularity at a — 0. One can find a limiting solution as a — > 0, in the following form 



E 



(m/£>) ^n!(n-l)!' 



(34) 



However we will get a real general solution in a compact form after integrating eq. (33) 
in the range — oo < k < +oo 



W{a,a,t) = 



\J AmDt 



(a) 1 exp 



ADt 



2l 



(35) 



In terms of p, it becomes 
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W(a,p,t) = 



; CXp 



2A-nGDpt 



(36) 



y/VWDGfPdH 

The two variable probability distribution function is Gaussian in cr, and diverges as 
a — > 0, where classical approach fails and quantum theory takes over. Now we write 
the expression for the transition probability 



P(a,a,t | a',a',t') 



^JAtiD (t - 1') 



cxp 



/M2l 



[(a - a') - (t - t')] 
AD (t - 1') 



(37) 



In terms of p and a it becomes 



cxp 



P (a, p, t | a, p, t') = (aa 



2AnGpp'D(t-t') 



(38) 



yjAnD (t - t') 

which represents the transition probability for the variables to change from the state 
(a',p') to (a,p) . Thus the scale factor a together with the density p evolves in a non- 
deterministic way, which in turn strongly influence the formation of large scale structure 
in the universe, since the evolution of the density perturbations are also depend on w. 
In all these cases we get Gaussian distributions, which is sharply peaked initially but 
spread out with time. 



5. Conclusion 

In the preceding sections we have described a stochastic approach to cosmology as a 
modification to the deterministic evolution of the universe in the standard model. In 
section 2, we have shown that fluctuations in the creation rates are physical processes 
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which can lead to a stochastic equation of state. A fluctuating w factor, in turn, 
will lead to fluctuations in the time - evolution of the energy density of the universe, 
as well as in the expansion factor for the universe. Thus both parameters become 
stochastic quantities, instead of remaining deterministic variables. A fluctuating w 
factor will also lead to fluctuations in the time-evolution of the Hubble parameter; i.e., 
the expansion rate of the universe becomes a stochastic quantity, instead of being a 
deterministic variable. We predict that such a stochastic expansion rate might have led 
to a randomness in the recession velocities of objects and will produce a scatter in the 
Hubble diagram. Using the redshift-magnitude data /x for the Type la supernovae of 
Perlmutter et al [4], which corresponds to their Fit C, and the technique used in Riess 
et al [5], we have computed the probability distribution function p(H \/i ) for a flat 
universe and compared it with the theoretical distribution [7], showing that both curves 
agree very well with D = 3.77 x 10 13 s and a half width of 0.011. A similar estimate 
is possible also with the distribution function of the density parameter, provided we 
are given a known data for the density parameter. Also if we have some explicit 
examples of models where a stochastic w emerges, the predicted value of D may be 
compared with our estimation, but in this paper, we have not made any attempts in 
this regard. Thus, if w is a fluctuating quantity, the evolution of the early universe 
becomes stochastic or non-deterministic (the scatter in the Hubble diagram indicates 
this), and the dynamical equations of those epochs are Langevin type equations, where 
one can evaluate the probability distribution functions of the variables. In [23], Berera 
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and Fang describe dynamically how, the stochastic fluctuations arising from various 
dissipations, generate seeds of density perturbations in the early universe, apart from 
the quantum fluctuations of the standard inflationary model. Here we attempted to 
make a stochastic approach to the early universe due to fluctuations in the mean 
equation of state, which is a classical phenomenon. Since the stochastic equation of 
state, lead to fluctuations in the time - evolution of the total density of the universe, 
the density contrast will also be fluctuating. We propose to undertake this study in 
detail in a future publication. 

To conclude, we note that the stochastic approach presented above is a modifi- 
cation to the standard model, when fluctuations are present. We have formulated a 
stochastic model and developed a set of non-deterministic, Langevin equations for the 
cosmo logical parameters in those epochs, under the assumption that the factor w(t) 
is a fluctuating quantity, when the universe is approximated by a many component 
fluid. It is expected that diffusion coefficient is a crucial factor in the evolution of 
the universe, especially in the early phase, where it influences the time - evolution of 
density parameter and scale factor for the universe. As fluctuations die out with time 
(D — > 0), the evolution becomes deterministic. 
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